Accuracy Control and Performance Enhancement of Linear Solvers for the Integrated Water Flow Model
نویسندگان
چکیده
In this paper, we describe the accuracy control and performance enhancement of linear solvers for the Integrated Water Flow Model (IWFM). This model is used by the State of California Department of Water Resources to assess the impact of climate change on water resources and the analysis of different conjunctive use scenarios across California. IWFM simulates groundwater, surface water and surface-groundwater interaction using an implicitly formulated Galerkin finite element approximation of the groundwater head in a multi-layer aquifer system. The computational efficiency of the simulation is governed by the efficiency of linear solvers for sequences of large-scale sparse linearized systems of equations. We firstly understand how multi-layer aquifer flow and stream-groundwater interaction affects the scaling, conditioning and sparsity structure of the linear systems. These properties guide the choice of scaling which, together with preconditioning, not only offset the ill-conditioning effects of multi-scale flow, but significanly improve the control of the linear solver forward error. Improved error control ensures that the accuracy of the solver is consistent with the accuracy of the initial data. We implemented a preconditioned Krylov subspace linear solver based on the Generalized Minimum RESidual (GMRES) algorithm and incomplete LU preconditioners and demonstrate how scaling improves forward error control in IWFM. We also performance benchmarked the new linear solver against the SOR method, a classical stationary iterative linear solver used in IWFM, and find an overall 7.7x speedup for the largest tested dataset. Further performance profiling shows that the new linear solver removes a major performance bottleneck in IWFM for the other datasets. Department of Computer Science, Kemper Hall, UC Davis, CA 95616 Corresponding author email: [email protected] Bay-Delta Office, Department of Water Resources, 1416 9th Street, Sacramento, CA 95814
منابع مشابه
Adaptive Accuracy Control of Nonlinear Newton-krylov Methods for Multiscale Integrated Hydrologic Models
In the popular Newton-Krylov methods for solving large-scale systems of nonlinear equations, inner linear systems resulting from outer Newton linearization are solved by Krylov iterative linear solvers. The accuracy control of Krylov solvers are based on the progress of the Newton iteration to achieve good local convergence while avoiding over-solving. In practice, the efficiency and robustness...
متن کاملViscous Models Comparison in Water Impact of Twin 2D Falling Wedges Simulation by Different Numerical Solvers
In this paper, symmetric water entry of twin wedges is investigated for deadrise angle of 30 degree. Three numerical simulation of a symmetric impact, considering rigid body dynamic equations of motion in two-phase flow is presented. The two-phase flow around the wedges is solved by Finite Element based on Finite Volume method (FEM-FVM) which is used in conjunction with Volume of Fluid (VOF) sc...
متن کاملSelection of Intermodal Conductivity Averaging Scheme for Unsaturated Flow in Homogeneous Media
The nonlinear solvers in numerical solution of water flow in variably saturated soils are prone to convergence difficulties. Many aspects can give rise to such difficulties, like very dry initial conditions, a steep pressure gradient and great variation of hydraulic conductivity occur across the wetting front during the infiltration of water. So, the averaging method applied to compute hydraul...
متن کاملA New Two Dimensional Model for Pollutant Transport in Ajichai River
Accurate prediction of pollution control and environmental protection need a good understanding of pollutant dynamics. Numerical model techniques are important apparatus in this research area. So a 2500 line FORTRAN 95 version code was conducted in which using approximate Riemann solver, couples the shallow water and pollution transport agents in two dimensions by the aid of unstructured meshes...
متن کاملEnhanced Predictions of Tides and Surges through Data Assimilation (TECHNICAL NOTE)
The regional waters in Singapore Strait are characterized by complex hydrodynamic phenomena as a result of the combined effect of three large water bodies viz. the South China Sea, the Andaman Sea, and the Java Sea. This leads to anomalies in water levels and generates residual currents. Numerical hydrodynamic models are generally used for predicting water levels in the ocean and seas. But thei...
متن کامل